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Abstract 

We perform the isospin decomposition of proton and neutron SLAC data 
for moderately low x in the region 0.01 < x < 0.1. The isovector part is 
well described by a power behaviour x a , where a does not correspond to the 
intercept of single a\ Regge trajectory, as expected. However, the observed 
power leads to the validity of Bjorken sum rule and it is consistent with the 
power extracted from all previous data using NLO evolution. At the same 
time, the isoscalar part behaviour may be interpreted as a partial cancellation 
between a positive non-singlet contribution and a strongly negative singlet 
one. Further experimental consequences are mentioned. 



I. INTRODUCTION 

The very accurate measurement of the neutron spin-dependent structure function g™ [|TJ , 
whose results have been presented recently, possess two remarkable properties. First, the 
values of neutron structure function at Q 2 = 5GeV 2 are rather large and negative in the 
region of moderately low x. Second, the data can be rather accurately fitted by the power 
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function x~ ' 8 and this power seems to be unexpectedly large. It is not obvious, where this 
large number, for which there is no indication in the proton data, is coming from, so it may 
affect the extrapolation to x = and cast some doubts on the validity of Bjorken sum rule. 

In the present paper we perform the isospin decomposition of the data for both proton and 
neutron. As a result, we conclude that the isovector contribution is well approximated by the 
power behaviour found earlier by an elaborate method based on NLO evolution |2| . It may be 
interpreted either as the manifestation of ln 2 x terms |§ or as the contribution of the Regge 
cut, produced by the a\ meson and a BFKL pomeron which has a high intercept. At the same 
time, for the isoscalar part one has the signature for a rather singular singlet contribution, 
compatible with the similar behaviour predicted by QCD [[|, although statistical errors are 
still rather large. This may imply a significant gluon polarization in the nucleon in this 
range of x, which could also clearly show up in double helicity asymmetries All at RHIC, 
if they turn out to be larger than estimated earlier p. 

II. BJORKEN SUM RULE VALIDITY 

The main feature of the new neutron data flj is large and negative g™ ~ — gf, measured 
with a good accuracy up to small x, say x ~ 0.01. The SLAC proton data || at Q 2 = 3GeV 2 
are positive and of roughly the same magnitude, but on the contrary, rather flat in this 
region. This could be understood, qualitatively, as a result of the interplay of a negative 
contribution at low x, responsible for the singular behaviour of the neutron, and a positive 
contribution at larger x. To check this assumption it is instructive to consider the isovector 
contribution. Since there is no clear evidence of scaling violations between Q 2 = 3GeV 2 
and Q 2 = 5GeV 2 in any polarized deep inelastic scattering experiment, we may neglect, for 
the time being, the effects of QCD evolution, because we are only interested in the results, 
provided by the data at the present level of statistical accuracy. 

It is then possible, by combining the SLAC neutron data with the SLAC proton data, 
to determine the quantity g^~ n = g\ — g™, entering the Bjorken sum rule, with a higher 
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accuracy, than for proton alone. This is because the neutron data are fortunately negative, 
so the difference is larger in magnitude than <jf, while the errors are practically the same for 
g\~ n and g±, given the better accuracy of the neutron data. 

Performing this analysis of the data, we see immediately a behaviour less sharp than 
that of neutron. Since phase space effects, as powers of (1 — x), are not so important in the 
region under consideration, we were looking for a simple power parametrization of the data, 
implied by Regge pole behaviour, namely 

gr n (x)=gr n (x )(-Y a , (1) 

\Xq J 

and we found that this works rather well, for 0.016 < x < 0.125 with the following choice of 
parameters: 

gT n {x) = O.U7x~ OA \ (2) 

as shown in Fig.l. The power we obtain is significantly smaller than the expected contribu- 
tion of the ai(1260) meson trajectory (~ 0.14). As a result, the contribution to the Bjorken 
integral from the region < x < 0.125 is large 

/■0.125 

/ dxg{- n {x) = 0.085. (3) 
Jo 

The region of higher x corresponds to a neutron contribution much smaller than that of 
the proton, the latter also providing a large contribution to the integral, and we have 

f 1 dxg{- n {x) « f 1 dxg{{x) = 0.09. (4) 

Vo.125 JO. 125 

The total contribution to the Bjorken sum rule, 

f 1 dxg{- n (x) = 0.175 (5) 
Jo 

appears to be in the fair agreement with the theoretical value. Note that when the final 
neutron data will be available, as well as more precise proton data, it will allow one a more 
serious analysis, taking also into account the effects of QCD evolution. 
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At the present moment we would like to stress, that by combining the current neutron 
and proton data, one is led to good agreement with Bjorken sum rule. Let us stress that 
this power is compatible with the one obtained using next-to-leading order (NLO) fit |2j 
to all previous data (—0.56 ± 0.21). It is consistent with the result of ln 2 x summation ||. 
We may also understand the observed power a = —0.45 by considering the intercept of the 
Regge cut associated to the a\ meson and a Pomeron, namely 

a = 1 — ap — a ai (6) 

provided we use the famous BFKL Pomeron |7j] with the intercept ap ~ 1.6 and a ai = —0.14. 

III. ISOSCALAR CHANNEL AND THE SINGLET CONTRIBUTION 



Since we found, that the sharp neutron structure function is not seen in the difference 
between proton and neutron, it should be attributed to the isoscalar channel. Also, the 
partial cancellation, we suspect to be at the origin of a flat proton structure function, 
should be manifested in this channel as well. To check this, we calculated the quantity 
gP+n ^ gP _|_ ^ re£ jiy g^Q-ws a ra ther flat structure for x > 0.035. The relative errors are 
much larger in this case because g\ +n is small due to the fact that g\ and have opposite 
signs. Of course this fact also implies a small value of the deuteron structure function in 
this kinematic region, which is barely consistent with the existing data [Q. 

This flat structure should be related to the interplay of the negative sharp contribution, 
showing itself in the fit x -0 8 to the neutron data [[|], and a positive contribution with 
a smaller power, dominating at larger x. Since there is no counterpart for such a sharp 
behaviour in the conventional Regge analysis, we make a strong, but natural assumption. 
Namely, we suggest that it is manifested in the 5'[/(3)-singlet channel. 

It is in this channel that a strong mixing between polarized quarks and gluons provides 
the anomalous gluon contribution to the first moment of g\ JTD|. Recent studies show that 
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generalized anomalous gluon contribution appears also in all moments [p]T|,[12| . At low x the 
quark-gluon mixing provides a strong correction to the subleading behaviour |Q, producing 
a power close to 1. One might expect that a similar effect is also present in the non- 
perturbative region where it can give rise to a strong x dependence in gluon distribution at 
low Q 2 , which is an initial condition in the approach just mentioned above. An example 



of such a non-perturbative contribution is given by instanton effects |13j which, however, do 
not provide yet a reliable quantitative estimate. 

Moreover, the SU(3)- nonsinglet part receives the contributions from the /i(1285)- and 
r/(547)-mesons trajectories. The first one has an intercept close to that of a\ occurring in 
the isovector channel, while the second one produces a smoother behaviour like a; 0,3 . For the 
first estimate we neglect the latter and find that the data are well described by the formula: 

g P+ n ( x ) = 0.145x-°- 45 - 0.03x-°- 87 , (7) 

as seen in Fig. 2. This formula is suggesting that the isoscalar contribution is approximately 
equal to the isovector one, which is not so surprising, in order to have the neutron structure 
function, dominated by the most singular power only. This would mean, that in this region 
of x, say between 0.01 and 0.1, one has 

Au(x) - Ad(x) ~ Ad(x) - As(x), (8) 

requiring a strong negative s-quark polarization. Apart from the difficulties of incorporating 
this result to current models of nucleon structure, it could also conflict with the Bjorken 
sum rule for the decay of strange baryons, implying that the integrals of both sides of 
this equation should be of opposite signs. Although the suggested above equality may be 
valid only in a limited region of x, and violations of SU(3) symmetry may be possible, it 
is more likely, that the r) contribution makes the x dependence of isovector and isoscalar 
combinations, different in the region under consideration, so we will have 

9i + ( x ) = C ns x 0,45 + C v x ' 3 — C s x a " . (9) 
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The numerical analysis shows, that the present inaccurate data, mainly for x < 0.035 
are equally well described with a wide range of coefficients Cf and C v . Bearing in mind the 
problem with the second Bjorken sum rule mentioned above, one may suspect a negative 
value for C v in order to reduce the isoscalar integral at larger values of x. On the other 
hand, the parameters of the singlet contribution are rather stable, C s ~ 0.03, a s ~ 1, which 
is again related to the good accuracy of neutron data. 

Note that, if this neutron behaviour is really an isoscalar phenomenon, as suggested 
by our analysis, one should observe a decrease and sign change for the proton structure 
function not too far from x ~ 0.005 at low Q 2 . This is the first major check of this result. 
However, negative proton structure function for much smaller values of x and low Q 2 have 
been considered in the literature and a sign change may also come at large Q 2 from the 
effect of QCD evolution [|. 

Note that the obtained power is also compatible with the NLO fit result, but the 
accuracy of the neutron data would allow to reduce the error. 

It is, of course, too early to relate unambiguously the observed behaviour to the results 
of ||, because of the limited experimental accuracy and some theoretical problems. In 
particular, it is not absolutely clear, to what values of x and Q 2 the results of || should be 
applicable. Nethertheless, the relative closeness of the experimental and theoretical numbers 
may be a signal, that the low x asymptotic behaviour is manifested rather early, especially 
in the neutron case, where it is not screened by a large nonsinglet contribution, like in the 
proton case. Note that clear evidence for negative no longer requires the negative gluon 
polarization, as was guessed in [|J, relying on earlier data. Moreover, we found that the 
formula (3.24) of || is not incompatible with the data, if some mean values (like suggested 
in [|J) of parton distributions are taken as an input. However, this approach does not allow 
one to extract the x-dependent parton distribution, and we shall use now the continuity 
with the region of average x in order to get an estimate. 

It is not clear, to which extent one should attribute a small x singlet contribution to 
quarks or to gluons. However, according to @] the contribution of gluons is dominant, 
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so we neglect the quark contribution in the present approach. Requiring the qualitative 
continuity in transition between low and average x, and applying the parton-like formula 
for the anomalous gluon contribution it seems natural to expect that the gluon distribution 
is behaving like 

AG(x) ~ x' - 87 . (10) 

Note, that this formula is assuming the simple relation between singlet contribution to g\ 
and AG at average x 

g{{x) = -^AG(x), (11) 



which is, strictly speaking, is valid for the first moment only. More generally, one has [14 



»«*> = -(£/. -r^y*- (12) 

where E(y) is the coefficient function, describing the gluon-photon interaction. Recent 
studies, based on the non-local generalization of the axial anomaly [IT ], support the following 
choice (1T2 



E(z) = 2(1 -z), (13) 

leading, for the gluon distribution of the type const x x~ a , to the relation 

2 n 

9{{x) = ~——-^AG(x). (14) 
a(a + 1) 07r 

As a result, for a = 0.87, AG(x) should be multiplied, for a given g{{x), by a factor 
~ 0.8 and the integral of AG(x) in the range 0.01 < x < 0.1 is about 1. 

The double helicity asymmetry All i n prompt photon production in pp collisions is 
directly related to AG(x) and for the center of mass energy a/s = 500GeV, which will 
be reached at RHIC, one is probing precisely this kinematic region of x. So given such a 
strong gluon polarization, we anticipate a larger A L l than previously predicted 0. Similar 
comments can be made for All in inclusive jet production. 



IV. DISCUSSION AND CONCLUSIONS 



We have performed here an analysis based on the present level of the experimental 
accuracy, still not enough to see the effects of QCD evolution, but providing interesting 
information about isospin structure. For this reason we restrict ourself to the SLAC data, 
and neglect the effects of QCD evolution in our analysis. Note that the SMC data jOJ are 
unfortunately not accurate enough to be used for such a simple isospin decomposition. The 
results on gf are certainly not incompatible with our g\ +n but in order to obtain g\~ n , it 
is necessary to extract from gf after subtracting g{, which enhances substantially the 
statistical errors. However, the elaborate statistical analysis using NLO evolution lead to 
the very similar powers, although the error for the singlet case is still very large. 

The presented simple picture of the nucleon structure is based on the two observations. 

i) As suggested by the E154 Collaboration, the g™ behaviour is well described by ~ a: -0 ' 8 . 

ii) From our simultaneous analysis of proton and neutron data, there is no indication of 
such a behaviour in gi~ n . Instead, it is well described by the ~ x~°' 45 , which leads to a good 
saturation of Bjorken sum rule. 

Consequently, the existence of a strong negative isoscalar contribution is implied by these 
two facts. It seems rather well established, and leads to predict a negative g\ for x below 
0.005. 

Both the interpretation of the nonsinglet behaviour as a ln 2 x terms (or cut produced by 
the BFKL pomeron), as well as the relation of the sharp behaviour of the singlet contribution 
at low x, and even further, to a strong gluon polarization, can be considered more speculative. 
It would be an unusual coincidence, that two rather different aspects of small- x physics 
manifest themselves in the same physical quantity. However, these assumptions seem to us 
possible, and they will be either supported or disproved by future more accurate data which 
will allow to elaborate a better analysis of the problem. 

We are indebted to C. Bourrely, A.V. Efremov, J. Ellis, S. Forte and E. Hughes for 
stimulating discussions and valuable comments. O.T. is grateful to Centre de Physique 
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FIGURES 

FIG. 1. Comparison for g\~ n between the curve given by eq. (2) and the SLAC data refs. HH 
For g\ at x = 0.0165 due to the absence of SLAC data we used the SMC data ref. [15]. 



FIG. 2. Same as Fig.l for g\ n with the curve given by eq. (7) 
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